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A new time domain formulation for 
broadband noise predictions 

J. Casper* and F. Farassatt 

NASA Langley Research Center, Hampton. VA 


ABSTRACT 

A new analytic result in acoustics called “Formulation IB,” proposed by Farassat, is used to 
compute the loading noise from an unsteady surface pressure distribution on a thin airfoil in the 
time domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with 
the loading source term. The formulation contains a far field surface integral that depends on the 
time derivative and the surface gradient of the pressure on the airfoil, as well as a contour integral 
on the boundary of the airfoil surface. As a first test case, the new formulation is used to compute 
the noise radiated from a flat plate, moving through a sinusoidal gust of constant frequency. The 
unsteady surface pressure for this test case is analytically specied from a result based on linear 
airfoil theory. This test case is used to examine the velocity scaling properties of Formulation IB 
and to demonstrate its equivalence to Formulation 1 A of Farassat. The new acoustic formulation, 
again with an analytic surface pressure, is then used to predict broadband noise radiated from an 
airfoil immersed in homogeneous, isotropic turbulence. The results are compared with 
experimental data previously reported by Paterson and Amiet. Good agreement between 
predictions and measurements is obtained. Finally, an alternative form of Formulation IB is 
described for statistical analysis of broadband noise. 


NOMENCLATURE 

A - correlation area (m 2 ) 

b = airfoil semi-span (m) 

C = airfoil chord (m) 

c 0 = ambient sound speed (m/s) 

/ = frequency (Hz) 

/ = geometry function for airfoil surface (Figure 1) 

E* = complex combination of Fresnel integrals (Equation (9g)) 
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Research Center, Hampton, VA 2368 1 

^Senior Research Scientist, Aeroacoustics Branch, Mail Stop 461, NASA Langley Research Center, 
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velocity-to-pressure transfer function 

time independent kernel function (Equation (24)) 

(o/U, convective wave number ( m 1 ) 

AC/2, reduced frequency 
streamwise integral length scale (m) 
correlation length (m) 

— y 

V/c 0 , Mach number vector 

— y — y 

M ■ rlr, Mach number in radiation direction 

— y 

M ■ v, Mach number in direction of v 
unsteady airfoil surface pressure jump (Pa) 
unsteady airfoil surface pressure (Pa) 
sound pressure radiated to observer (Pa) 

surface pressure gradient in the direction of surface velocity (Pa/m) 
length variable associated with radiation distance in retarded time 
(Equations (23b) and (23c)) 
correlation function of w (Equation (19g)) 

— y — y 

x — y, sound radiation vector (m) 

power spectral density of w in A-space (m 3 /rad-s 2 ) 

observer time (s) 

uniform freestream speed (m/s) 

unsteady streamwise velocity (m/s) 

airfoil velocity vector in reference frame fixed to undisturbed medium (m/s) 
unsteady upwash velocity (m/s) 

[x p x 2 , x 3 ] r , observer position in reference frame fixed to undisturbed 
medium (Figure 1), fixed to airfoil reference frame in Section 5 only 
[y , , y 0 , 0] 7 , surface source position in reference frame fixed to undisturbed 
medium (Figure 1), fixed to airfoil reference frame in Section 5 only 
VI — M 1 

cJf acoustic wavelength (m) 

Mk/p 2 

directivity angle (Figure 7) 

unit inward facing normal on surface edge (Figure 1) 

angle between airfoil surface normal and radiation direction (Figure 1) 
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p () = ambient density (kg/m 3 ) 

t = source time (s) 

<I> = random phase variable 

u) = circular frequency (rad/s) 

Subscripts 

1, 2, 3= chordwise, span wise, and normal coordinate directions (Figure 3) 

1A = solution calculated by Formulation 1A 

IB = solution calculated by Formulation IB 

ret = evaluated at retarded time t — r/c Q 

1. INTRODUCTION 

The development of analytical methods to predict noise radiated from an air flow over 
a rigid body has been an active research topic within the aeroacoustic community for 
decades. Research in this area has, in large part, been motivated by the desire to 
incorporate the results of aeroacoustic analysis into an aerodynamic design 
methodology. The present work is similarly motivated, and the resulting formulation 
should lend itself well to an engineering design tool suite. 

The current work is specifically focused on the calculation of far field noise that 
results from fluctuating pressure on a solid surface. The acoustic analogy [1] provides 
a framework for developing methods to predict noise from many types of sources, 
including noise due to unsteady surface loading. Such noise is mathematically 
described by the loading source term, or “dipole term,” of the Ffowcs Williams- 
Hawkings (FW-H) equation [2], Because the noise due to an air flow over a rigid 
surface is typically dominated by dipole radiation, the acoustic formulations of interest 
in this work are determined by solutions of the FW-H equation with the loading source 
term, i.e. neglecting the thickness and volumetric terms. 

The solution of the FW-H equation can be written in many forms. A new solution, 
called “Formulation IB,” is presented herein. Most types of broadband noise, including 
noise that is dominated by an airfoil’s leading and trailing edges, can be calculated with 
the proposed Formulation IB, which is the simplest loading noise prediction formula 
known to date. This simplicity makes the new formulation highly suitable for statistical 
analysis of broadband noise for rotating surfaces. 

In Section 2, Formulation IB is derived for the specific case of a flat plate 
undergoing a general nonuniform motion. For low Mach numbers and distant 
observers, the dominant term in this formulation is a far field surface integral that 
depends on the time derivative and the gradient of the airfoil surface pressure. The 
formulation also contains a contour integral on the boundary of the airfoil surface. This 
line integral vanishes along the trailing edge if the Kutta condition is imposed. 

In Section 3, Formulation IB is used to calculate the noise radiated from a flat plate 
moving through a sinusoidal gust of constant frequency. The unsteady surface pressure 
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used in this test case is an analytical result from linearized airfoil theory that is taken 
from the work of Amiet [3, 4]. A mesh refinement study is performed to demonstrate 
the equivalence of Formulation IB with Formulation 1A [5], a previously developed 
acoustic formulation that is also a solution of the FW-H equation. Results from this 
analytical test case are also used to examine the velocity scaling properties of 
Formulation IB, which are found to be consistent with the results of Curie [6] and 
Ffowcs Williams and Hall [7]. The directivity of the noise induced by a periodic gust 
is also examined. 

In Section 4, the single-frequency surface pressure in Section 3 is extended by 
spectral representation to serve as an analytic broadband source model for incident 
turbulence noise. This stochastically modeled surface pressure is used in Formulation 
IB to predict broadband noise caused by the interaction of a NACA 0012 airfoil and 
homogeneous turbulence. The resulting calculations are compared with experimental 
data previously reported by Paterson and Amiet [8]. In Section 5, an alternative 
acoustic formulation is described for statistical analysis of broadband noise. 

2. ACOUSTIC FORMULATION 

— > 

Consider a flat, finite surface moving in the plane x, = 0 along a velocity vector V . The 
velocity vector V and the surface geometry are related to the coordinate axes as 
pictured in Figure 1. Let f(x v x 2 , t) denote a geometric function where / = 0 on the 
surface edge and/ > 0 on the interior of the surface. Let v = denote the unit normal 
which lies in the plane of the surface, is normal to the edge, and is directed inward 
(Figure 1). Note that V need not be constant in space or time. The only stipulation on 
the velocity is that the motion of the surface is in the same plane as the surface. 



Figure 1 . Schematic for the derivation of Formulation IB. 
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Let x = [Xj, x 2 , x 3 ] 7 denote the position of an observer, and y = [_y 3 , y 2 , 0] 7 denote 
the position of a source point on the plate’s surface (Figure 1). The spatial frame of 
reference is fixed to the undisturbed medium. The unsteady surface pressure p(y, t) 
gives rise to sound that initiates at source time t and radiates to the observer along r = 
x* — y*. This sound is described by p'(x, f), the perturbation pressure that arrives at the 
point (x v x 2 , x 3 ) at time t. For some of the subtle mathematical details that are involved 
in the following derivation, see [9], 

From the FW-H equation, the loading noise is given by a solution of 

i ay 

— -t* - vy = - ^ • [ pn H(f ) S(x 3 )] [la] 

c o clf 

where c Q is the ambient sound speed and/i is the unit surface normal that, for the present 
case, is equivalent toe 3 , the unit vector in the direction of the x 3 -axis. The function H(f) 
is the Heaviside step function and S(x 3 ) is the Dirac delta function. Evaluating the 
divergence in Equation (la) yields 

i ay 

— - vy = -p(x v X 2 , t) H(f) S'(x 3 ) [lb] 

c o dt 

where S'(x 3 ) denotes differentiation with respect to x 3 . 

Equation (lb) is the wave equation with an inhomogeneous source term, and its 
formal solution in an unbounded domain requires the Green’s function 8(g)/4Trr, where 
g - t — t + r/c 0 , and t and rare the observer and source times, respectively. The solution 
of Equation (lb) can then be written in the form 

/•! f 8(g) 

4 73 y (X ,t)= - J J^ 3 ~^p(y v y 2 , T) H(f) 8'(y 3 ) df dr 


Let t — » g and integrate with respect to g. The result can be written 


4 7 rp'Ot, t) = -J^ y (p] ret H(F) S'(y 3 ) df [2a] 

where the subscript “ret” denotes evaluation at retarded time t — r/c Q , and F is 

F(y v y 2 , x , t) =f(y v y 2 , t - r/c Q ) = [/] ret [2b] 


Integration with respect to y 3 on the right-hand side of Equation (2a) yields 
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[3a] 


Performing the differentiation in the integrand of Equation (3a) yields 
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dy 3 
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ip] 


ret 


H(F) 


— [f 3 p] ra H(F) + — 
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3 /\ 


[3b] 


r 3 P 


ret 


8(F) 


where p is the time derivative of pressure evaluated by an observer that is fixed with 
respect to the medium at rest, and r 3 is the third component of the unit radiation vector 
r - ~r!r. Clearly, then, r, = e 3 ■ r = cos 6, where 0 is the angle subtended by the surface 
normal and the radiation vector, as shown in Figure 1. The first and second terms on 
the right-hand side of Equation (3b) are of the similar form Q FI ( F ). In [10], it is shown 
that the integration of these two terms in Equation (3a) can be written 


L Q II (F) dy x dy 2 = f Qdl = f Q dS [4a] 

■'ll? 2 J F>0 v >0 [1 - M.] ret 

where rfX is the surface area element of the acoustic planform of F > 0 . Also, /W = ,vl • r 

is the Mach number in the radiation direction, where A? = [7/c 0 is the local Mach 

number vector of the surface. 

The integrated value of the third term on the right-hand side of Equation (3b) is 
determined as follows. This integral is of the form 

J=[ q(y t ,y2)8(F)dy t d yi [4b] 

The differential surface element dy x dy 2 can be written [10] 

dSJlF 

dy l dy 2 = d£jBi= ^ [4c] 


where c/fAtnd t/iA^are differential elements of arclength that are, respectively, parallel 
and normal to the the surface edge defined by F = 0 as shown in Figure 2. The notation 
denotes the surface gradient in the y t vyplane. Moreover, it can be shown [10] that 
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dC d * 
W\ = [1 - M ,Jret 


where M is an element of arclength along the surface edge defined by / = 0. Equation 
(4b) can then be written 


r _ dlJF c qyy v y 2 ) * 

= L q(y v y 2 ) S(F > ^ & = I. ^ djj= I 

% 2 j f = o J f= o 


q(y v y 2 ) 


q(y v y 2 ) 


\v 2 f\ 


|v/l 


J f= 0 [1 - Ml 


[4e] 


H ret 


Note that the surface time derivative d/Ydr in Equation (3b), and contained in q(y v y 2 ) 
in Equation (4e), is referenced to the undisturbed medium. However, df/dr can be 
related to the local normal in-plane velocity of the edge contour as 


Df 3/ ^ 3/ , 

= — +V-Vf= — +V-v=0 

Dt , „ 3t 3r 


If =o 


[5a] 


where the subscript “/ = 0” denotes evaluation on the edge contour. This relation yields 
df /3t = — P ■ v. 


Figure 2. 


(F > 0 ) 







Differential surface element in Equation (4c). 


All three terms on the right-hand side of Equation (3b) are now integrated in 
Equation (3a) over the physical surface/ (x v x 2 , t) > 0 by using Equations (4a) and (4e). 
Note that in Equation (la), p is referenced to the medium at rest, e.g., as measured by 
a transducer that remains stationary as the surface passes by it. The quantity p can also 
be related to dp/dr, the time derivative of pressure in the reference frame of the moving 
surface, e.g., as measured by a transducer attached to the surface. 
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This relation is 


dp dp 

dr ds 


[5b] 


where dp/ds is the gradient of p in the direction of P, V is the local magnitude of P, 
and ^ is in the direction of the velocity P of the surface in the reference frame fixed to 
the undisturbed medium. 

Incorporating all of the above results into Equations (3a) and (3b), the solution of 
Equation (la) can be written 


4tt //(x\ t) = [ 

Jf- 

L 


7>o 


p cos 6 


( dp/dr — V dp/ds) cos 6 
c n r (1 - Ml 


dS 


ret 


r 2 (1 - Ml 


r 

M v p cos 0 

dS - 


O 

II 

r( 1 -M) 


di f. 


ret 


[ 6 ] 


where M = M ■ v, the Mach number in the direction of v. Equation (6) will be referred 
to as Formulation IB. 

The first and third integrals in Equation (6) are far field terms, and the second 
integral is the near field term. The relative contributions of the terms in Equation (6) 
are now considered, under the conditions of low Mach number and an observer in the 
acoustic far field, i.e., 


M«l, t»li [7] 

where X is a typical acoustic wavelength of interest. With respect to M and r, the 
surface far field integral, i.e., the first integral in Equation (6), is proportional to 1/r, 
whereas the second and third integrals are proportional to 1 hr and Mir , respectively. 
Therefore, the far field surface integral dominates the signal under the conditions in 
Equation (7). 

Again, Equation (6) is valid for the case of nonuniform surface motion. Therefore, 
Equation (6) can be used, as is, to predict loading noise from rotating surfaces. Its 
predecessor. Formulation 1A [5], is signicantly more complicated in its rotational form 
and cannot be approximated by only one surface integral in the far field. Such a 
signicant simplication for far field calculations makes Formulation IB more suitable 
for statistical analysis of broadband sources for rotating surfaces. A statistical 
formulation based on Equation (6) will be addressed in Section 5. However, the focus 
of the current work is the application of Formulation IB in the form of Equation (6), 
as will be demonstrated in the following two sections. 
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3. SINUSOIDAL GUST OF CONSTANT FREQUENCY 

Any noise prediction made with Equation (6) will be only as good as the input surface 
pressure p(y, r). Such time -dependent pressure data could be taken from experimental 
measurements or a computational fluid dynamics (CFD) calculation. However, in this 
initial work, an analytic expression is used for p(y , r) to serve as input data for the new 
acoustic formulation. To this end, an analytic formulation from thin airfoil theory is 
used to describe the unsteady surface pressure that results from a sinusoidal gust of 
constant frequency. This particular surface pressure formulation is chosen as an 
example to establish the equivalence of Formulation IB and Formulation 1A [5] and to 
examine the velocity scaling properties of Formulation IB. In addition, the frequency 
dependence of the far field directivity is discussed. 

3.1. Surface pressure from thin airfoil theory 

Consider a rectangular flat plate, in rectilinear motion, as in Figure 3. The velocity 
vector V = [ — C/, 0, 0] T , where U is a constant subsonic speed. For the following 
examples, the plate’s surface and its boundary, / > 0, are defined by the rectangle {0 < 
Xj < C} X { —b < x 2 < b] in the plane x 3 = 0. This surface is assumed to have an 
unsteady pressure distribution that is analytically prescribed from thin airfoil theory, as 
discussed below. 



< c > 

Figure 3. Schematic for the constant-frequency loading noise problem in Section 
3; an infinitely thin rectangluar surface in rectilinear motion. 

Amiet [3, 4] presents closed form expressions for the unsteady pressure on the 
surface of a thin airfoil of infinite span. The airfoil is assumed to move rectilinearly 
through a sinusoidal gust. Analytical methods are used to solve the two-dimensional, 
time-dependent linear potential equation by representing the solution as a product of 
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spatial and temporal solutions. The solution is represented as a truncated series in 
which higher order terms are neglected. 

A complex-valued representation for the airfoil surface pressure is assumed to arise 
from a stationary gust in one spatial dimension. This gust can be written in the 
stationary variable x } — Ut as 

w(x j - Ut) = w Q e- ,k(x < - u,) [8] 


where k = u>/U is the streamwise convective wave number, and vv Q is the gust amplitude. 
This gust and the airfoil surface pressure that it generates are, for now, considered 
functions of a single amplitude and frequency. 

The unsteady pressure that arises at a point (y p y 2 , 0) on the airfoil surface as a result 
of the gust in Equation (8) can be written 


APCy,, 0 = 2p Q Uw Q g(y v k ) e ikut 


[9a] 


where p Q is the ambient density, k - kCH is the reduced frequency (based on the semi- 
chord), and g(y , , k) is a transfer function whose form is dependent on the frequency of 
interest. The factor of two in Equation (9a) indicates that the pressure is assumed to be 
antisymmetric between the upper and lower surfaces, and this expression thereby 
accounts for the pressure on both sides of the airfoil, i.e., the pressure jump. In [8], the 
suggested parameter to del ineate between the low and high frequency regimes is p = 
Mk/(3 2 , where = I — M 2 , with a demarcation value of p = 0.4. 

For low frequencies, p <0.4, the transfer function is 


8(y v V= — 


c 

— - 1 

y\ 


1 

2 


G s (k*)e'**«b „ , p< 0.4 


[9b] 


where k* -k/f3 2 , G s is the classical Sears function [11], which, for the present work, is 
approximated by 


G s (k*) = 


1 

1 + 2.4k* 


+ 277k* 


[9c] 


as suggested in [8], and 

q(y v M) = M\2y x IC - 1) + (1 - /3) In M + j8 ln(l + f3) - In 2 [9d] 

For high frequencies, p > 0.4, the transfer function is the sum of a leading edge 
solution and a trailing edge correction [4], i.e., 


g(y,, k) = (gj + g 2 ) e 


*[2/a(1 - M)y\!C + it! 4] , ^ > q 4 


[9e] 
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where 


1 


c?i(y P k) = 


[2tt^j( 1 + M)/C]i 


[9f] 



-1 + (1 + i) E* [4/r( 1 - y,/ Q] 



and 


E*(Q = ‘ T du - m - iS(£) 

J o (2 ttm)J 


[9g] 


The quantities C(£) and S( £) are the Fresnel cosine and sine integrals and will be 
numerically evaluated by the formulas in [12]. The final representation for the unsteady 
surface pressure p(y v t), assumed to be a real quantity, is 


The pressure jump is negative in Equation (9h) because the acoustic formulation 
was derived from Equation (la), in which the unit surface normal e 3 is assumed to point 
into the fluid, i.e., in the positive x, direction on the upper surface, and in the negative 
x 3 direction on the lower surface. Therefore, using the same positive surface normal on 
both sides of the airfoil, the sum of the pressure on both sides is A P = F* U pp er ~ ^lower 
and this expression is the negative of the conventional notion of a pressure jump. 

Some differences are noted here between Equations (9e) and (9f) and their 
counterparts in [3] and [8]. Such differences include the choice of the coordinate-axes 
origin and the spatial normalization employed by the authors. Another difference arises 
from the use of complex conjugates in the present paper, which serves to make the 
present notation more consistent with the eventual broadband formulation in [8], which 
is employed by the present authors in Section 4. The above surface pressure was 
originally proposed in order to derive an expression for unsteady lift that was ultimately 
incorporated into a frequency domain acoustic formulation [3, 8, 13]. However, in the 
present work, the unsteady pressure jump itself will be used as input to Equation (6) for 
a time domain prediction. Finally, note that for both transfer functions, the pressure jump 
vanishes at the trailing edge x = C; i.e., the Kutta condition is satised. The correctness 
of this trailing edge boundary condition has been a subject of debate for many years, for 
example as reviewed by Howe in [14]. However, as the response of the airfoil leading 
edge is expected to dominate the noise that arises from an incident gust, determining the 
most appropriate treatment at the trailing edge is not an issue for the current problem. 


p(y v T) = ^{-AP(yj, t)} 


[9h] 


3.2. Grid refinement study 

The surface pressure in Equations (9a) - (9h) is now used to numerically demonstrate 


218 


A new time domain formulation for broadband noise predictions 


the equivalence of Formulation IB and and the loading noise terms of Formulation 1A 
[5]. Formulation 1A forms the basis of WOPWOP, a rotor noise prediction code 
developed at NASA Langley Research Center [15]. To predict loading noise from an 
airfoil in uniform rectilinear motion. Formulation 1 A is simplied to 


The entire formulation here is integrated on the surface interior. At first glance, the 
form of Equation (10) appears no more complex than Equation (6), but only because 
of the simplicity of uniform, rectilinear motion. Unlike Equation (6), Equation (10) is 
not applicable to a rotating surface. The loading noise terms of the full Formulation 1 A 
are indeed applicable to rotational flow, but the full formulation is more complex than 
Equation (10) (see [15]). 

The far field noise radiated from a thin airfoil in a one-dimensional, single- 
frequency sinusoidal gust is now calculated with Equations (6) and (10). Let p' lB and 
p\ A denote the sound calculated by Equations (6) and (10), respectively. If the input 
surface pressure p(y, r) is exactly specied at every point on the airfoil surface, then the 
only error in the numerical solution of Equations (6) and (10) arises from the 
quadrature formula that is chosen to perform the surface and contour integrations. In 
this case, the midpoint rule is the quadrature of choice. Therefore, the equivalence of 
Equations (6) and (10) is demonstrated if the difference \p' w — p' ] A \ diminishes in mesh 
refinement like the cumulative error expected from the midpoint rule, i.e. that the error 
is proportional to the square of the mesh spacing. 

The plate’s rectangular dimensions are determined by a chord length of C = 0.5 m 
and a span of 2b = 2.0 m. The plate is moving at a Mach number of 0.2, and the sound 
speed is taken to be 343 m/s. The ambient density is p () = 1.23 kg/m 3 , and the upwash 
amplitude is w () = 0.05 U. The observer position for this test case is x = [—1,0, I ] 7 , in 
meters. Figure 3 roughly depicts the direction of this observer position, although the 
distance is not to scale. 

The calculation is performed for one time period of the surface pressure uctuation 
at frequencies of/= 25 Hz and /= 1 kHz, with 32 time steps in each period. These 
choices of frequency, at the prescribed observer location, will test both the near field 
and far field equivalence of the two formulations. Note that the transfer functions in 
Equations (9b) and (9f) are singular at yq = 0, and the spatial derivative of Equation (9f) 
is singular at y x = C. Although both singularities are integrable in an analytic sense, 
they would cause the quadrature error to deviate from that of the midpoint rule which, 
by its denition, requires sufficient smoothness throughout the interval of integration. 
Therefore, the domain of streamwise integration is restricted to an interval of the form 



[ 10 ] 
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eC <y l < (l — e) C [11] 

where £ is a small, positive parameter. 

Each calculation is performed on a sequence of six surface grids: { 10 X 40], [20 X 
80}, {40 X 160}, [80 X 320}, {160 X 640}, and {320 X 1280}. Because the surface 
pressure is cast in only one spatial variable y v and the observer location is symmetric 
relative to the airfoil span, the acoustic prediction is relatively insensitive to the 
discretization in y T These discretizations are employed as a more general example of 
the types of resolution that would be required for a case in which spanwise pressure 
variation would signicantly affect the far field acoustics. 

U 



Figure 4. Second level grid (20 X 80) for convergence study in Section 3.2. 



Figure 5. Grid refinement validation for equivalence of Formulations 1 A and IB. 
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The grid is clustered via Roberts’ transformation [16] near the leading and trailing 
edges of the plate in order to accommodate the parameter e = 0.02 on the coarser 
meshes. The 20 X 80 grid is shown in Figure 4. The maximum values of | p' lB — p’ ]A \ 
during each time period are shown as a function of the number of grid points on a log- 
log plot in Figure 5. The abscissa N c is the number of surface elements in the chordwise 
direction. The dashed line represents a fictitious quantity whose values are specifically 
calculated to be directly proportional to N ( , 2 . The slopes of both calculations are 
visibly parallel to a slope of —2, thereby demonstrating that 


which is the cumulative error expected from the midpoint rule of integration. 

3.3. Velocity scaling laws 

Attention is now turned to the way in which the intensity of the far field noise, as determined 
by Equation (6), will scale as a function of velocity when the surface pressure is described 
by Equations (9a)-(9h). The conditions of a far field observer in a low Mach number flow 
(Equation (7)) will be assumed. The Mach number range of interest is 0.01 < M < 0.2. 
Velocity scaling laws will be determined for both compact and noncompact sources. 

First, the source is assumed to be compact, i.e., C <s A. This condition will be 
achieved by requiring the frequency to be proportional to velocity, / ~ U/C, for a 
sufficiently low range of frequency. This assumption is consistent with the compact 
source argument put forth by Curie [6], If a frequency range of 4 Hz < / < 80 Hz is 
chosen for the Mach number range 0.01 < M < 0.2, then the acoustic wavelength A 
ranges from 85.75 m at M - 0.01 to 4.29 m at M - 0.2, for an ambient sound speed of 
343 m/s. If the plate’s physical dimensions are as in the above mesh refinement 
problem (C = 0.5 m), then the ratio A/C ranges from 8.58 at/= 80 Hz to 171.5 at /= 4 
Hz. These ratios are assumed sufficient for an acoustic source to be considered 
compact. Note that p < 0.4 throughout these ranges of Mach number and frequency. 

The calculation is performed on a 100 X 400 surface grid, with the streamwise 
integration interval restricted as in Equation (11), with sufficient grid clustering near the 
leading and trailing edges to allow for £ = 0.003. The observer position is chosen at a 
distance of 100 m, directly above the plate’s center, i.e., x = [0.25, 0, 1 00] 7 i n meters. 
This location places the observer in the acoustic far field for the entire range of frequency. 

The upwash amplitude is w () = 0.05 U for each of 50 equally spaced Mach numbers 
between 0.01 and 0.2. A separate calculation is run for each Mach number and its 
corresponding frequency. The surface pressure in Equations (9a)-(9d) and (9h) is used 
as input to Equation (6) to predict the far field sound p'{x , t). Each calculation is 
performed for one acoustic period T of the corresponding frequency, with 64 time- 
steps. The average intensity I(x) of the acoustic signal at the observer ?, assuming 
spherical spreading, is then calculated by 


\p\b - P\A = °Wc 2 ) = 0{ M) 
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The average acoustic intensities for a compact source, as a function of Mach 
number, are represented as squares in Figure 6. The slope of these results on a log-log 
plot can be visually determined by observing their proximity to the dashed line with the 
slope of six. This U 6 proportionality is consistent with Curie’s result [6], as expected 
from the conditions placed on the calculations. 

The demonstration of a velocity scaling law is next desired for a noncompact source 
where /i > 0.4. Therefore, the restriction that / ~ U/C must be lifted so that / is 
independent of U. The simplest such condition is that /is constant, in which case the 
Strouhal number remains a function of U. A series of calculations is again performed, 
as above, with the only parameter change being that the frequency is held constant at 
1 kHz throughout the range of Mach number. The input surface pressure for this case 
is given by Equations (9a) and (9e)-(9h). The computed acoustic intensities at 100 m 
are represented by circles in Figure 6. In this case, the acoustic intensity scales 
approximately as J/ 5 , a result that is consistent with Ffowcs Williams and Hall [7], 



Figure 6. Velocity scaling properties of Formulation IB. 

3.4. Directivity 

As a final exercise in this section, the directivity of a single -frequency source is 
examined. The radiated noise p'{x , t) is calculated at many locations on a circular arc in 
the plane x, = 0 that is geometrically centered on the plate’s upper surface, as shown in 
Figure 7. The arc trajectory (r, iff) is determined by r = 3 m and 0 < i//< it. The directivity 
is determined by the peak pressure amplitude |p'| x that is calculated at each position 
on the circular arc during one period in time for a given frequency. The flat plate’s 
dimensions and surface discretization are the same as in the previous example, and there 
are 128 time-steps in a period. The observer path, 0 < ijj < tt, is discretized into 128 
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equally spaced locations. The freestream Mach number is 0.2 and the gust amplitude is 
w () = 0.05 U. Figure 8 shows the results, in polar form, for frequencies of 1 kHz and 2 
kHz. The higher frequency results are scaled to show both loci on the same plot. As 
expected, the 2-kHz solution has twice as many lobes as the 1-kHz solution. 
Furthermore, the general downstream directivity of the lobes is consistent with the 
results of previous authors e.g., [17, 18]. Finally, note the frequency dependent positions 
of the lobes with respect to a fixed observer. For example, an observer at <// = 90° will 
receive a signal that is near the peak amplitude for the central lobe of the 1-kHz signal, 
whereas the 2-kHz signal is near a local minimum for the same observer. This frequency 
dependent characteristic of directivity is discussed further in the next section. 




\ 

\ 



2 

Figure 7. Schematic for directivity calculation; observer on circular path in plane x 2 = 0. 


90° 


/ = 1 kHz 
/ = 2 kHz 



I P Imax, 

Figure 8. Directivity for two sources of constant frequency; 2-kHz locus scaled by 
a factor of 2.5. 
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A. BROADBAND PREDICTION WITH COMPARISON TO EXPERIMENT 

The analytic surface pressure in the previous section is extended to model a broadband 
source on a slender airfoil at zero angle of attack. This broadband surface pressure is 
used as input to Formulation IB to predict far field radiation, and the results are 
compared with experimental measurements. Amiet [19] has previously proposed a 
broadband solution to this problem in the time domain by using the transfer functions 
in Equations (9e)-(9g). His resulting solution was a Fourier transform of the frequency 
domain solution in [8] and [13]. In the current work, the high frequency formulation in 
Section 3 is explicitly extended to a broadband source application with user-specied 
spectral content. 

A.1. Experiment description 

The experiment that is modeled in this section is reported by Paterson and Amiet [8], 
A NACA 0012 airfoil is placed between two vertical plates at zero angle of attack in 
the test section of an open jet wind tunnel. A schematic of this experimental setup is 
shown in Figure 9. The airfoil has a chord length of 0.23 m and a span of 0.53 m. 
Turbulence is generated by a grid upstream of the airfoil. Noise propagates from the 
test section into an anechoic chamber that is instrumented with six microphones. The 
microphones are located on the tunnel centerline, on an arc of radius 2.25 m, relative 
to the airfoil’s geometric center. The microphone location of interest here is at 90 
degrees (Figure 9). Far field noise measurements of the incident turbulence on the 
airfoil are determined by subtracting microphone measurements, with and without the 
model in the test section, at each of five tunnel speeds: 40, 60, 90, 120, and 165 m/s. 


Figure 9. 



Schematic for incident turbulence experiment in [8]. 


A.2. Broadband analysis 

For prediction purposes, the airfoil is modeled as a flat plate in a periodic gust which 
gives rise to an unsteady surface pressure that is a broadband extension of the analytic 
formulation in Section 3. The airfoil geometry is oriented with respect to the coordinate 
axes as in Figure 3, with [0 < x x < C] X { —b < x 2 < b}, where C = 0.23 m and 2b = 
0.53 m. The components k v k v and of the wave number vector are oriented along 


224 


A new time domain formulation for broadband noise predictions 


the axes x v x 2 , and x y respectively. As encountered by the airfoil surface in the plane 
x., = 0, the normal component of the turbulent velocity field can be written 


where w(k v k 2 ) is the wave number spectrum of the gust amplitude, defined as the 
inverse Fourier transform of w(x v x 2 , t). Furthermore, this two-component velocity 
spectrum is interpreted as a fully three-dimensional spectrum whose form results from 
the integration of all wave numbers that are normal to the airfoil surface: 


The complex-valued, unsteady surface pressure jump that arises at a point (y p y n , 0) 
on the airfoil surface from the incidence of a turbulent velocity field of the form in 
Equation (12) is given by 


A straightforward approach for predicting the desired broadband far field 
measurements is to use the real part of Equation (13) as input to Formulation IB. The 
velocity spectrum required in this approach can be derived from well known power 
spectrum formulas for homogeneous, isotropic turbulence, such as von Karman’s 
formula [20]. A formulation for the dual wave number transfer function g(y p k v k~) is 
given by Adamczyk [21], 

However, because one of the objectives of the current work is to reproduce the 
results of Paterson and Amiet [8] from a time domain perspective, an approach similar 
to that taken in [13] will be taken to model the surface pressure. Amiet [13] argues that, 
within certain limitations, integration over all spanwise wave numbers is not required. 
His conclusion, derived mathematically in the frequency domain, is that only one 
spanwise wave number contributes to the sound detected by an observer in a given 
location. In particular, Amiet focuses on an observer in a spanwise symmetric location, 
for which only the zero spanwise wave number needs to be considered. This result is 
argued to be exact in the limit of infinite span and a good approximation for an airfoil 
of finite span that encounters a high frequency disturbance. 

The derivation of Amiet’s analytical result can be generally described as follows. First, 
Equation (13) is transformed into Fourier space, with the stationary variable x { — Ut 
being replaced by to = k l TJ. Then, a two-point cross-correlation function is formed and 
related to the far field power spectrum through KirchhofPs formula [22] and Curie’s 



[ 12 ] 
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result [6]. In order to follow a similar line of reasoning in the time domain. Equation 
(13) itself must be related to the far field acoustic pressure through Formulation IB. In 
the case of a distant observer directly overhead of a finite-span airfoil, the terms r, M ., 
and 0 in Equation (6) are weak functions of y x and y 2 , and therefore will be considered 
constants. Furthermore, for the observer position considered here, the differences in 
retarded time, as a function of airfoil surface location, can be neglected. These 
assumptions are consistent with the acoustic model employed by Amiet [13]. For the 
present problem, including the above assumptions. Equation (6) is approximated by 


(x, t ) = 


cos 0 

c 0 r(l -M.) 



d 

—p(y v y v t - r/c 0 ) + 


U — p(y v y 2 , t - r/c 0 ) 
^1 


M cos 0 pi 

dy 2 dy \ + — I P(°> yv 1 ~ r /c o)dy 2 

r(l — M ) J -b 


[14] 


where the over-bars ond, r, and M. denote mean values over the airfoil surface, and 
therefore the retarded time t — r/c 0 is constant for fixed t. The near field term (second 
integral in Equation (6)) has been deleted because of its 1/r 2 proportionality. However, 
the third integral must be retained because the Mach number is not necessarily assumed 
to be small. Only the leading edge line integral survives the contour integration in the 
third term in Equation (14) because M = 0 along the airfoil tips and the Kutta condition 
is assumed along the trailing edge. 

Before the surface pressure p is specied. Equation (14) is further simplied. For 
convenience, the terms 6 and M . will be neglected, as they are small (M ~ 0 and cos 0 
~ 1) for a distant observer directly above the airfoil. With these additional 
simplications, if —A P in Equation (13) is substituted for p(y v y 2 , t — r/c Q ) in Equation 
(14), the far field acoustic pressure can be approximated in the form 


2 pJJ rc r b pc*. 

47 rp'{x,t)~ — J J J J (F, (y ! , k ! , k n ) e l(k < Ut k ^^dk 2 dk , dy 2 dy , 

[15a] 

2p,,U M p b poo poo 

+ _ !L J b J J J 2 (0, k v k 2 ) e l,k < ut - * 2 * 2 )dk 2 dk l dy 2 


where 


T^ y v k v k 2 ) = -w(k v k,) i k { U gOq, k k 2 ) + U — g(y v k v k 2 ) 

dy i 


g — i r/c 0 


J 2 (0, k v k 2 ) = -W(k v k 2 ) g{ 0, k v k 2 ) e~ ik t f/c o 


[15b] 
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Sufficient conditions [23] on A P and its derivatives have been assumed for the 
commutation of integration and differentiation in Equations (15a) and (15b). The y 2 
integration in Equation (15a) can now be explicitly evaluated, yielding 


47777' (x , t) ■■ 


2 p a U pc 


in: 


2sin(k 2 b) 


c 0 r ° k 2 


7 x iy v k v k^e 1 k ' Ut dk 2 dk ] dy l 


[16] 


2 P 0 U M roo poo 2sin(k 2 b) 


a 


J 2 (0, k v k 2 )e' k "‘dk 2 dk 1 


Integrating with respect to k 2 , the term sin (k^b)/k 2 acts like a Dirac delta function 
when integrating over an unbounded domain, and the result is 


47777' (x , t) ■ 


2p,U 


2pJJ M 


J J 2ir f 1 (y v k v O)e ik i Ut dk l dy l 

J 2ttJ 2 ([), k v Q)e ik ' u, dk ] 


[17] 


Equation (17) indicates that only the zero span wise wave number contributes to the 
noise detected by the far field observer. Equations (16) and (17) are time domain 
analogies to Equations (15) and (17) in [13]. Furthermore, Equation (17) suggests that, 
as input to Equation (6), the surface pressure jump in Equation (13) can be simplied to 


A P(y v t ) : 


2irp Q U 


w(k v 0)g(y p k v 0 )e'^ u, dk 


[18] 


Note that the y 2 integration will be explicitly performed when Equation (18) is input 
to Equation (6). Only the k , integration is neglected. Furthermore, the evaluation of the 
surface pressure terms at retarded time t — r!c {] will be executed in an exact fashion, as 
prescribed by Equation (6). Note also that, as a function of observer time t, the surface 
pressure in Equation (18) is constant in —b < y 2 < b. Although the present analysis 
involves an airfoil of finite span for which end effects are expected, Amiet [12] argues 
that these end effects are negligible except at very low frequencies. Therefore, the 
surface pressure in Equation (18) will be applied to the entire airfoil surface, including 
the regions near the tips. 

The evaluation of the surface pressure in Equation (18) is accomplished by first 
recognizing the turbulent fluctuations as a stochastic process. This process can be 
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approximated by a truncated series whose coefficients are chosen such that the 
autocorrelation of the series forms a Fourier transform pair with its power spectrum 
(see, e.g., [24]). This property is achieved by evaluating the spectral amplitudes w(k v k n ) 
as a function of the power spectral density (PSD) of w(k v k 2 , t). To this end, the infinite 
wave number domain, — °° < k l < °°, in Equation (18) is integrally discretized and 
truncated such that k { _ N < k ln < k l N . The largest wave number k { N represents an 
“upper cutoff’ wave number, beyond which the spectral density amplitude is 
considered negligible or is out of range of experimental measurement. The unsteady 
surface pressure jump in Equation (18) is then approximated by 

N 

A P(y v t ) « 2np 0 U^ 0 e*. g(y v k l n , 0) e‘ k ^< 

n=—N 

k ln = n A k v n = 0, ±1, ±2, ..., ±N [19a] 

M i = k J N 


The phase angles { } are independent random variables uniformly distributed on 
[0, 27 t]. The gust amplitudes {A 0 ] are related to the PSD of the turbulence and will be 
discussed momentarily. The transfer function g(y p k { , 0) is the limiting case of 
Adamczyk’s [21] dual wave number transfer function and is given by 

g(y v k ln ,0) = (g l + g 2 )e- i ^nO-M)y l i c + * / 4-k l J [19b] 

where g l and g 2 , now functions of k ] , are otherwise identical to Equation (9f), 


SlOb’ Kr) = 


1 

[2Trk Xnyi {\+M)lcf- 


s 2 (y v k \J = 


-1 +(1 +0 £*[4/x„(l - y,/C)] 
[27Tk u {\+M)i 


[19c] 


and E* is the same complex combination of Fresnel integrals as in Equation (9g). The 
low frequency transfer function is not used in these broadband predictions because the 
experimental facility is anechoic for frequencies above 200 Hz, and the parameter p, is 
greater than 0.4 at this frequency or above for all five tunnel speeds. 

The gust amplitudes [A n0 ] are evaluated by 

A = [.S' ( k . , 0) M,]s 

«,0 L vpw v !,«’ y 1 J 


[19d] 
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where S (k v k 2 ) is interpreted as a three-dimensional PSD of w, into which all normal 
wave numbers are integrated, 



S (k., k~, kk) dk-. 

wvtA 1’ 2’ y 3 


Although A k 2 does not appear in Equation (19d), the expression in Equation (19a) 
still represents a dimensionally consistent surface pressure with respect to its role as 
input to the acoustic formulation in Equation (6). This fact arises from the previous 
acoustic analysis in which k 2 was explicitly integrated in order to obtain the expression 
for the input surface pressure jump in Equation (18). Moreover, the acoustic analysis 
in Equations ( 14)-( 17) alleviates the need for the scale factor that was previously 
reported in [25]. It is important to note that the broadband predictions to be discussed 
in the following subsection are generated without any reliance on a scale factor or other 
arbitrary constant. 

As this grid generated turbulence is assumed to be homogeneous and isotropic, the 
PSD can be related to von Karman’s formula [20]. The formula for S mv (k v k 7 ) used in 
the present calculations is derived from von Karman’s energy spectrum in [13]. 


5 (k.,kj= 

U>VtA 1 ’ // ,7 


4 u 2 k 2 +k 2 


9l T k 2 [1 +k 2 +k 2 ]i 


where u 2 is the streamwise mean-square turbulence, and 


k t Jn T(D 

k. = — , i = 1,2, k 

1 K e Ci m 


[19e] 


[19f] 


and is the streamwise integral length scale 



R (x,)dx, 

wwA V 1 


where R w is the upwash correlation function defined by 


1 |* T 

R ww (r) = {w(x v t ), w(x j + r, t)j = lim — J w(x p t) w(x 1 + r, t) dt 


[19g] 


Note that the denition in Equation (19g) implies that the stochastic process is 
assumed ergodic. 


aeroacoustics volume I ■ number 3 ■ 2002 


229 


The two-component formula in Equation (19e) was derived in [13] by integrating 
the von Karman energy spectrum over all T, components. The k 2 component is then set 
to zero for the present calculations. Values for Tr and 2l| are determined by 
experimental measurement. Fink [26] reports that the turbulence intensity that results 
from the grid used in this experiment can be approximated by the empirical formula 




U 


0.04 


u 


-0.2 


[19h] 


where the reference speed is U ref = 60 m/s. Fink also reports a measured value for the 
integral length scale to be ^ = 3.2 cm. 

All of the above expressions and measurements are incorporated into Equation 
(19a). The final representation for the unsteady broadband pressure on the airfoil 
surface is then given by the real part of — A P in Equation (19a). This broadband surface 
pressure is used as input to Formulation IB to predict the acoustic pressure p'(x , 1) at 
the microphone location. 

Using symmetry arguments and algebraic manipulation, the indicial bounds for the 
surface pressure’s spectral representation are altered so that the domain includes only 
positive wave numbers. As input to Equation (6), the resulting real-valued surface 
pressure can be written 


p(y v T) = -4tt - p 0 t/j£ A n 0 [B n cos (k ln U T+ $ - a) + D n sin (k l n U t+ - a)] [20a] 


n= 1 


where the upwash amplitudes A Q are evaluated by Equation (19d)-(19h), and the 
quantities B n , D n , and a n are given by 


B „ = S i(y v k u) + 


C(i n ) + S(i n ) - 1 

[2irk u {\ + M)]5 


- ao 

[ 2 t 7* 1jB (1 + M)f 


[20b] 


Ti _ tt 
a = 2/i ( 1 — M ) — — k, + — 

n v 7 C 1 J 1 L 

where k t ) is evaluated in Equation (19c), and C(£ n ) and S(^ n ) are the Fresnel 
cosine and sine integrals in Equation (9g) with f = 4 /r ; ( 1 — y/C). The summation in 
Equation (20a) begins at n = 1 because S ( 0, 0) = 0, by Equation (19e). 
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M3. Time domain predictions 

The lower frequency bound, and therefore the fundamental frequency, for all five 
calculations is chosen to be 10 Hz. The upper frequency for the predictions is chosen 
according to the upper frequency for which measurements are available for each tunnel 
speed. For U = 40, 60, and 90 m/s, the upper bound is / v = 2.5 kHz. The upper bounds 
for U = 120 and 165 m/s are f N = 3.5 kHz and 4.5 kHz, respectively. For all five 
calculations, the numerical bandwidth is A/= 10 Hz. Each calculation is performed for 
one period of the lowest frequency, T - 0.1 s.The numerical solution is sampled at the 
Nyquist frequency, i.e., Ar = TUN, where N is the upper cut-off index in Equation 
(20a). The calculation is performed on a 100 X 230 surface grid with the streamwise 
integration interval restricted as in Equation (11), and with sufficient grid clustering 
near the leading and trailing edges to allow for £ = 0.003. As in the constant frequency 
case, because the surface pressure is cast in only the one spatial variable y v and the 
observer location is symmetric relative to the airfoil span, the acoustic predictions were 
found to be relatively insensitive to the discretization in y T Again, the current 
discretization is employed as an example for a more general case, such as when the 
spanwise differences in retarded time become important for an observer whose location 
is asymmetric relative to the airfoil geometry. 

The experimental microphone position for which comparisons are made is at a 
distance of 2.25 m from the model, and at an angle of 90 degrees relative to the 
geometric center of the model. The coordinate system for the calculation is such that 
the x 2 -axis coincides with the center span line, so that the microphone position is in the 
plane x 0 = 0,as in Figure 7. The measured observer position for the prediction is, then, 
x = [0.115, 0, 2.25] r in meters. 

The position of the microphone relative to the airfoil is corrected for refraction 
caused by the presence of a shear layer that forms downstream of the upper lip of the 
square nozzle exit and is positioned between the model and the microphone (Figure 9). 
This correction is based on geometrical acoustics with an assumption of a zero 
thickness shear layer, and is reported by Amiet [27]. Shear layer corrections that are 
based on such simple formulations are reasonable for the present case with the 
microphone directly above the source. The required correction in the microphone 
position is signicant. At a measured angle of 90 degrees, the corrected angles ranged 
from approximately 84.5 degrees for U - 40 m/s to 68.2 degrees for U = 165 m/s. 

In addition, the amplitude of the radiated noise is also corrected for the presence of 
the shear layer; however at a measured angle of 90 degrees the amplitude correction is 
not signicant, especially for the lower tunnel speeds. The computed sound pressures 
p'(x, t) were corrected by factors ranging from approximately 0.997 for U = 40 m/s to 
0.942 for U = 165 m/s. The microphone position is corrected for the shear layer in a 
preprocessing step. After the far field noise is calculated at the corrected position, the 
results are then post-processed for amplitude correction. In this way, the corrected 
predictions can be compared with the experimental results “as measured.” 

Figure 10 shows the far field signal p'(x , t ) that is predicted by Formulation IB, for 
one fundamental period in time, at the experimental microphone location, for a tunnel 
speed of 165 m/s. Shear layer corrections for amplitude and directivity are included in 
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this plot. In order to compare with experimental measurements in [8], the time domain 
results for all five tunnel speeds are Fourier analyzed and the resulting spectral 
amplitudes {P n } N =l , are used to compute the sound pressure level (SPL) spectrum of 
the far field radiation. These SPLs are determined by 


SPL(f n ) = 20 log 



n = 1, 2, ..., N 


[ 21 ] 


where the reference pressure is / J r( = 20/xPa. The SPLs are converted to a 1 Hz 
bandwidth by reducing the values in Equation (21) by 10 log(A f). This narrowband 
conversion is consistent with the experimental SPLs which were measured at a 
bandwidth of 55.7 Hz and reduced by 17.5 dB [8]. 



time , S 

Figure 10. Predicted far field signal, U = 165 m/s; microphone at 90°, 2.25 m above 
airfoil center. 

The predicted and measured far field SPLs for the five tunnel speeds are shown in 
Figure 1 1 . The experimental measurements in Figure 1 1 were obtained by digitizing the 
measurements plotted in Figure 4 of Reference 8. The solid symbols represent those 
measurements for which the difference between the noise with and without the airfoil 
model was considered too small, and these measurements are therefore subject to 
greater uncertainty. The agreement with the measured data is very good, particularly in 
the range of 200 Hz to 1500 Hz, where the spectra are peaked. The relatively signicant 
error in the prediction at the lowest tunnel speed is consistent with the low speed 
prediction in [8]. This error results from the fact that, as the air speed decreases, the 
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gust wave lengths for the lower frequencies that dominate the noise will decrease and 
eventually become of the order of the airfoil thickness, thereby making the flat plate 
approximation increasingly invalid. 



Figure 11. Predicted and measured far field noise spectra; microphone at 90°, 2.25 m 
above airfoil center; experimental data reproduced from [8]; solid 
symbols denote low signal-to-noise ratio. 

The existence of multiple local extrema in a spectrum is often interpreted as the 
existence of multiple noise producing mechanisms; however, in the present case, the 
multiple extrema result from the changing placement of lobes, as a function of frequency, 
relative to a fixed observer. This phenomenon was described in Section 3.4 (Figure 8). A 
similar trend in local extrema is also apparent in the experimental measurements at the 
two highest tunnel speeds, where the signal-to-noise ratio was much higher. 

In Figure 12, a comparison is made between the current time domain prediction and 
the frequency domain prediction of Paterson and Amiet [8], for the tunnel speed U = 
165 m/s. The highest speed is chosen for the large frequency domain for which the 
spectrum is predicted. The frequency domain prediction in Figure 12 was obtained by 
digitizing 51 points on the appropriate locus in Figure 4 of Reference 8. As the length 
of the square symbols is equivalent to 0.5 dB on the ordinate axis, clearly, the two 
predictions are very close throughout the entire frequency range, and nearly identical 
in the peak range of 200 Hz to 1500 Hz. That the current time domain predictions are 
so similar to the predictions in [8] is not surprising, as the current predictions rely upon 
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the same unsteady surface pressure formulation as those in [8]. However, the 
decoupling of the aerodynamics and the acoustics that is apparent in Equation (6) 
makes Formulation IB amenable to any type of input available, whether analytical, 
experimental, or computational. 



Figure 12. Comparison of predicted far field noise spectra, U = 165 m/s. 

5. STATISTICAL FORMULATION 

Often, when aeroacoustic experiments are performed, surface pressure correlations are 
extremely useful in characterizing noise due to an air flow over a model. Under certain 
conditions. Formulation IB is readily transformed into an expression that aids in the 
statistical analysis of broadband noise. Such an expression for the autocorrelation of 
the far field acoustic pressure is developed in this section. 

A far field autocorrelation formulation is derived for the case of a thin airfoil moving 
rectilinearly at constant velocity. The assumption of uniform velocity simplies the 
algebra considerably, and the extension to general motion is straightforward. The 
following derivation follows the statistical analysis for jet noise presented by Morris and 
Farassat [28]. In the case of an observer in the acoustic far field, r » A, the near field 
term (the second integral in Equation (6)) can be ignored because of its 1/r 2 
proportionality. Also, for notational convenience, the spatial and temporal derivatives of 
pressure in the first integral are written as a single time derivative of pressure,/), when 
evaluated in a reference frame that is fixed relative to the medium at rest, as in Equation 
(5b). In this fixed reference frame, the far field acoustic pressure can be written 


477 p' (X*, /) = [ 
J f> 0 


p cos 9 


c 0 r( 1 - M) 


dS - 

ret 



M v p cos 9 
r( 1 -M) 


di f . 
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Equation (22) is now transformed to the reference frame fixed to the airfoil, i.e. the 
wind tunnel reference frame. In this frame of reference, the emmision distance r, M has 
the following form: 


1 

r ret = TT l M ( x i ~ ^i) + R] [ 23 a] 

H 

where 


R = [(*! - y,) 2 + P\x 2 - y 2 ) 2 + y3 2 Jcf] 1/2 


Furthermore, 


[r (1 - M ; .)] ra = R 


and 


x 3 

[cos 0] ret = 

r ret 


P 2 X 3 

m (x y - yJ + R 


[23b] 


[23c] 


[23d] 


Therefore, in the wind tunnel reference frame, a time independent function 'J((x , v ) 
is defined by 


cos 0 (3 2 x 3 

= -x = ‘Klx, y) [24] 

r(l— M ( )] ra R[M(x f - >>,) + R] 

Equation (24) is an important result because it shows that, in the laboratory 
reference frame, the only time dependent quantities in Equation (22) are p and />. Thus, 
Equation (22) can be written in a form that shows the dependence of p'(x, t ) on T\ ~y, 
and t explicitly: 


1 p 

4777?' (x , 0 = — J 3 C(x , y ) p(y , t - r ret /c Q )dy s 
Cft j >0 



M v (y ) X(*, y) p(y, t - r ret l C 0 ) dP 


[25] 


where dy s denotes an integral surface element dy { dy T Using the notation (-,-) to 
represent an ensemble average, as in Equation (19g), assuming ergodicity, the 
autocorrelation of the acoustic pressure is evaluated as follows: 
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16 7 r(p'(x, t ), p'(x , t + t)) = 

-T L L {p(y,t- r ret /c 0 ), p(y\ t + T- r'Jc Q ))dy s dy' s 

I- I- f - r ret /c o)’ p(y'’ t+f ~ r 'm /c o )Ws df t 26 ] 

c 0 J f > 0 J f = 0 

- — j M M V «' (p&’*~ r ret /c o)’ p&> t + T- r W c o 'iWs 

+ J^_ Q L» M v K XX {p(y, t - r ret /c 0 ), p(y’, t + f- r' r X 0 ))d< ? df 

All primed quantities on the right-hand side of Equation (26) indicate replacement 
of the variable y with y', e.g., X = X& y’)- 

Now, let j? and y' be related by y*' = y + rf, where rf - [rj v 17 -,, 0] r is a vector in the 
plane of the airfoil. Furthermore, it is assumed that |r^ <K r, e| , in which case r' e[ can be 
approximated by 


'ret 


; r ret + V ■ V 


v r ret 


where 


V r : 
v ret 


M + 




R 


P 2 


yi 


R 


, 0 


[27a] 


[27b] 


Two autocorrelation functions, F ■ and F , and a cross-correlation function F are 

pp pp pp 

therefore defined by 


F pp(y\ v,t) = ( p(y , t), p(y+v,r + f)) 

F pp (y I T) = (p(y . , t), p(y + ij, T + r)} [28] 

7^(y; ^ f) = (p(y , t ), p(y +r),T + f)} 

Using Equations (27) and (28), the ensemble-averaged quantities on the right-hand 
side of Equation (26) have the following meaning: 

{p(y,t- r ret /c 0 ), p(y\ 1+f- rJc Q - rf- V/ ret /c 0 )} - F pp (y- y,f-p- V/ ret /c 0 ) 
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{p(y,t- r ret /c 0 ), p(y\ t + f- r^/c 0 - rf- V/ ret /c 0 )) = F pp (y; v,f-rf- V/ ret /c 0 ) [29] 

{p(y,t- r ret /c 0 ), p(y’, t + f- r' ret /c Q - rf ■ V/ ret /c 0 )} - F pp (y- p, f - rf ■ V/ ret /c 0 ) 

{p(y, t - r ret /c 0 ), p(y\ t + f- r'Jc Q - rf- V/ ret /c 0 )} = F pp (y- - p, -f + rf ■ V/ ret /c 0 ) 
Equation (26) can thus be rewritten using the above correlation functions explicitly: 
16 TT^fp'ly, t), p'(x , t + t)) = 

— J- dy s ‘Mx, y) | KL(x, y +p) F pp (y; p,f-p- V y r ret /c Q )dff 
C 0 V >0 A c 

- — dy s ‘K&, y) M v (y + rf) Kl{x, y + rf) F pp (y; -rf, -f + rf - V/ ret /c 0 )af/ [30] 
J._ d£ M v (y ) Kit y) J Kit, y + rf) FJy\ p,f -rf- V r^c^dff 

C q v - 0 A c 

+ j d£ MJy) Kit t \ MJS + v) d<i(ty +V) F pp (y; v,f~rf- V/Jc 0 )dP 

Here, A , and L are correlation area and correlation length, respectively. Clearly, 
performing the integration in the four inner integrals on the right-hand side of Equation 
(30) will result in four functions of t, y , and f. Calling these function F { to f 4 in order 
of their appearance. Equation (30) can be rewritten 


1 f 

16 Triply, t ), p'(x , t + T))= — Kjt, y ) [F^x, y, f) - c Q F 2 (x, y, f)]dy s 

c o 


1 f 

J~_ M „(y) ‘K(y, y ) [E 3 (^, y, f) - CqF ' 4 (x, , y, f)]dF 

c o f~° 


s 

[31] 


The usefulness of the analysis presented above depends on the knowledge of two 
autocorrelation functions and one cross-correlation function, all of which are expected 
to be obtained from experimental measurement. If the surface pressure fluctuations are 
assumed to be stationary in time, then only the correlation function F (y , rf f) needs 
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to be measured, as the other two correlation functions can be derived from it [29]. It is 
expected that, of the three correlation functions, F^(y ; rf, t) will dominate the acoustic 
pressure, and therefore, F x (x , y , f) will be the only signicant contributor to the right- 
hand side of Equation (31). Equation (30) indicates that, at a location y on the airfoil, 
only the turbulent pressure fluctuations over a correlation area near that point will 
contribute to the sound at the observer location 3c. Furthermore, the time dependence of 
the correlation functions is displaced by rj • V r re| /c 0 because of the retarded time effect. 
This displacement is dependent upon airfoil Mach number and the radiation direction 
3? — y (the visual direction). Perhaps as important as its direct application to the 
prediction of far field noise from surface pressure measurements. Equation (31) can be 
used for qualitative analysis of broadband and trailing edge noise, which is the subject 
of ongoing research. 

CONCLUDING REMARKS 

The prediction of broadband noise from rotating machinery and airframes is currently 
the subject of intense research in aeroacoustics. The physics of broadband noise 
generation are well understood as the result of the pioneering research of Howe [14, 30, 
31], Amiet and coworkers [3, 4, 8, 13, 19, 32], and Brooks and coworkers [33, 34]. The 
previous work of these aeroacousticians, and many others, has clearly demonstrated 
that any successful broadband loading noise prediction requires an understanding of 
two physical processes: the character of the time-dependent surface pressure that 
provides the acoustic source, and the manner in which that source gives rise to an 
acoustic signal. 

The most prevalent types of unsteady surface pressure associated with broadband 
loading noise are those that arise from the interaction of the airfoil with a gust or the 
passage of eddies near the trailing edge. Obtaining this fluctuating pressure distribution 
analytically, numerically, or experimentally is itself a difficult problem. For this reason, 
past researchers have most often resorted to modeling the surface pressure, using 
guidance from experiments to aid in the development of these models. Today, high 
resolution surface pressure fluctuations can be obtained from turbulence simulations in 
realistic situations where the airfoil geometry and kinematics are accurately modeled. 
Therefore, the improvement of the acoustic radiation model becomes an important 
research topic. In the past, acoustic radiation models were most often developed for 
airfoils in uniform rectilinear motion. In addition, other restrictive assumptions, such 
as far field positioning of the observer, were often used to simplify the acoustic 
analysis. 

The present work develops a simple and general acoustic result in the time domain, 
based on the solution of the loading noise term of the Ffowcs Williams-Hawkings 
equation. This new solution, called Formulation IB, is, to date, the simplest analytical 
result for the prediction of loading noise and is suitable for statistical analysis of 
broadband noise for a surface in general motion. The new formulation has been 
validated with time domain calculations that predict the noise due to incident 
turbulence on a NACA 0012 airfoil. The time domain predictions are found to be in 
excellent agreement with the frequency domain predictions of Paterson and Amiet [8] 
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as well as with their experimental measurements over a signicant range of tunnel speed. 
These results are, to the authors’ knowledge, the first successful broadband predictions 
in the time domain. 

Because of its relative simplicity. Formulation IB lends itself to a straightforward 
extension that provides an autocorrelation of the acoustic pressure in terms of two 
autocorrelation functions and one cross-correlation function of the surface pressure. 
Both numerical and experimental surface pressure distributions can be used in this 
result to yield information on broadband noise statistics. 

The authors advocate the use of time domain methods in the prediction of 
broadband noise. Because of the decoupling of the aerodynamics from the acoustics, 
the chief advantage of time domain methods is their potential for direct use of time- 
dependent surface pressure statistics from experiments or computer simulations in the 
prediction of broadband noise. 
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